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On the basis of a first-principles, relativistic electronic structure theory of finite temperature 
metallic magnetism, we investigate the variation of magnetic anisotropy, K, with magnetisation, 
M, in metallic ferromagnets. We apply the theory to the high magnetic anisotropy material, Llo- 
ordered FePt, and find its uniaxial K consistent with a magnetic easy axis perpendicular to the 
Fe/Pt layering for all M and to be proportional to M 2 for a broad range of values of M. For small 
M, near the Curie temperature, the calculations pick out the easy axis for the onset of magnetic 
order. Our results are in good agreement with recent experimental measurements on this important 
magnetic material. 

PACS numbers: 75.30.Gw, 75.10.Lp, 71.15.Rf,75.50.Bb,75.50.Ss 



By accounting for relativistic effects such as spin-orbit 
coupling on electronic structure, recent 'first-principles' 
theoretical work has been successful in describing trends 
in the magnetocrystalline anisotropy (MCA) of magnetic 
materials. Ql 0. Q This is useful for the understanding 
of permanent magnetic properties, domain wall struc- 
ture, magnetic nanostructures etc. One aspect, however, 
which has received scant attention from such ab-initio 
theories is its temperature dependence. Modelling this 
phenomenon seems to rest largely on the seminal papers 
by Callen and Callen and others published nearly 40 
years ago which focussed on the temperature dependence 
associated with single ion magnetic anisotropies. In this 
letter we investigate how far this approach can be justi- 
fied for metallic magnets on the basis of a first-principles, 
material specific, parameter-free theory of how the mag- 
netocrystalline anisotropy depends on temperature. The 
proposed theory involves a detailed, relativistic descrip- 
tion of the electronic structure and hence it includes a 
complete desciption of the spin-orbit coupling. The ther- 
mally induced magnetic fluctuations are taken into ac- 
count by a relativistic generalisation of the, by now well- 
tried, disordered local moment (DLM) picture. 0, ||| 

The topic has recently received extra impetus from 
extensive experimental studies of magnetic films and 
nanostructures and their technological potential. For 
example, the recent fabrication of assemblies of smaller 
and smaller magnetic nanoparticles holds considerable 
promise for the design of ultra-high density magnetic 
data storage media. 7] But this is hampered by a parti- 
cle size limit set so that thermally driven demagnetisa- 
tion and loss of data is avoided over a reasonable stor- 
age period. This limit can be lowered by using mate- 
rials with high magnetocrystalline anisotropy, K, since 
the superparamagnetic diameter of a magnetic particle 
is proportional to (ksT/K)^ , where ksT is the ther- 



mal energy. |8( In this context, the chemically ordered 
LIq phase of equiatomic FePt, which has high uniaxial 
MCA (4-10 10 7 ergs/cm 3 or up to 1.76 meV per FePt 
pair has been subjected to much attention and 

arrays of FePt nanoparticles with diameters as small as 
3 nm have been synthesised. 0, EI A way to write to 
media of very high K material is by temporary heat- 
ing- EE 1 1 ■ J-l] The magnetocrystalline anisotropy is 
reduced significantly during the magnetic write process 
and the information is locked in as the material cools. 
Modelling this process and improving the design of high 
density magnetic recording media therefore requires an 
understanding of how K varies with temperature. 

So for the first application of our theory we have cho- 
sen Llo-ordered FePt. Given its technological potential 
there have recently been some careful experimental stud- 
ies of its fundamental magnetic properties. 0, flil Il5j | 
These show a strong temperature dependence to K. We 
find good agreement with these data. In particular we 
find K{T) oc (M(T)/M(0)) 2 over a broad magnetisa- 
tion range, in line with the results reported in two of 
these three studies. Jl2> llj| Notably the low tempera- 
ture behavior is qualitatively different from that of the 
single-ion anisotropy models used over many years. 

The magnetic anisotropy of a material can be conve- 
niently expressed as K = k 1 g 1 (n) where the fc 7 's are 
magnetic anisotropy constants, n is the direction of the 
magnetisation and <? 7 's are polynomials (spherical har- 
monics) of the angles 9, <f> describing the orientation of h 
and belong to the fully symmetric representation of the 
crystal point group. For a uniaxial fcrromagnet such as 
Llo-ordered FePt, K = k + fc 2 (cos 2 - 1/3) + ■ • •. 

As the temperature rises, K decreases rapidly. The 
key features of the results of the early theoretical work 
on this effect Q are revealed by a classical spin model 
pertinent to magnets with localised magnetic moments. 
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The anisotropic behavior of a set of localised 'spins' 
associated with ions sitting on crystalline sites, i, in 
the material, is given by a term in the hamiltonian 
H an — Yli X) 7 kjg 7 (si) with §i a unit vector denoting 
the spin direction on the site i. As the temperature 
is raised, the 'spins' sample the energy surface over a 
small angular range about the magnetisation direction 
and the anisotropy energy is given from the difference 
between averages taken for the magnetisation along the 
easy and hard directions. If the coefficients fc 7 are as- 
sumed to be rather insensitive to temperature, the dom- 
inant thermal variation of K for a ferromagnet is given 
by K(T)/K(0) = (gi(s)) T /( 9l (s)) The averages (• ■ -) T 
are taken such that (s)t = M(T), the magnetisation of 
the system at temperature T, and I is the order of the 
spherical harmonic describing the angular dependence 
of the local anisotropy i.e. I = 2 and 4 for uniaxial 
and cubic systems respectively. At low temperatures 
K(T)/K(0) « (M(T)/Af(0))^ +1 )/ 2 and near the Curie 
temperature T c , K[T)/K{Q) w (M(T)/M(0)) 1 . These 
features are borne out rather well in magnets where the 
magnetic moments are well-localised, e.g. rare-earth and 
oxide magnets, but it is questionable whether such an 
analysis should hold for itinerant ferromagnets. Q Here, 
we examine these issues in the context of the Llo-ordered 
FePt intermetallic compound for which careful experi- 
ments [HEi find K(T)/K(0) = (M(T)/M(0)) n ,where 
n = 2 instead of n = 3, over a large temperature range. 
As will be shown presently, our ab-initio calculations are 
in good agreement with this surprising result. 

Magnetocrystalline anisotropy is caused largely by 
spin-orbit coupling and receives an ab-initio description 
from the relativistic generalisation of spin density func- 
tional (SDF) theory. Q Up to now calculations of the 
anisotropy constants K have been suited to T = OK only. 
They treat the spin-orbit coupling effects using either 
perturbation theory 0] or a fully relativistic one @, ■ 
Typically the total energy, or the single-electron contri- 
bution to it (if the force theorem is used), is calculated 
for two or more magnetisation directions, hi and h 2 sep- 
arately, i.e. E ni , E n2 , and then the MCA is obtained 
from the difference, AE. Since the MCA is a small part 
of the total energy of the system, in many cases of the or- 
der of /J,eV, it is numerically more precise to calculate the 
difference directly [II] so that AE = - J Efi [N(e; hi) - 
N{e-h 2 )]de-\n{E F2 -h 2 ){E Fl -E F2 ) 2 + 0{E Fl — Ep 2 ) 3 
where Ep , Ep are the Fermi energies when the system is 
magnetised along the directions hi and h 2 respectively 
and n(s; h) and iV(e; h) are the density of states and 
integrated density of states respectively We have used 
this rationale with a fully relativistic theory to study the 
MCA of magnetically soft, compositionally disordered 
binary and ternary component alloys Ea. Il7j , the effect 
upon it of atomic short-range order [2| and recently its 
dependence in FePt upon the extent of Ll long range 
chemical order [lif . 

There are a number of calculated values of the MCA 



of completely Llo-ordered FePt at T = OK in the lit- 
erature ranging from 1.2 to 3.9 meV (7-22 10 7 ) 
ergs/cm 3 . K for FePt has been measured at room tem- 
perature to be 6.6 10 7 ergs/cm 3 although more recent 
measurements suggest that this value could exceed 
10 8 ergs/cm 3 for completely ordered FePt. The easy 
axis in both experiment and theory is along the c-axis, 
(0,0, 1), perpendicular to the Fe and Pt layers. 

Our calculations of the MCA of materials use spin- 
polarised, relativistic multiple scattering theory and an 
adaptive mesh algorithm 20] for Brillouin zone integra- 
tions such that the numerical precision of our calcula- 
tions is to within 0.1 /xeV. 0, 0] These attributes are 
also important for the theoretical calculations of the tem- 
perature dependence of the MCA described below. We 
calculate the MCA of ordered FePt at T = OK to be 
1.696 meV. We start from a self-consistent field, scalar 
relativistic calculation (atomic sphere approximation) of 
the electronic structure and effective potentials for the 
Fe and Pt sites. We then perform a further fully rela- 
tivistic electronic structure calculation, recalculate the 
Fermi energies Ep and Ep and determine the MCA. For 
Llo-ordered FePt we find that if Ep and Ep are artifi- 
cially changed by a small amount, the MCA varies sig- 
nificantly. For example, the MCA jumps from 1.696 to 
2.751 meV if the E F 's are lowered by just 0.2 eV. This 
sensitivity may explain, in part, the range of published 
values of the MCA of FePt. We can also deduce that 
the magnetic anisotropy of this strong permanent mag- 
net might be further enhanced by replacing a few atomic 
percent of Pt with Ir. 

In a metallic ferromagnet at T = OK the electronic 
band structure is spin-polarised. With increasing tem- 
perature, spin fluctuations are induced which eventually 
destroy the long-range magnetic order and hence the 
overall spin polarization of the system's electronic struc- 
ture. These collective electron modes interact as the 
temperature is raised and are dependent upon and affect 
the underlying electronic structure. For many materials 
the magnetic excitations can be modelled by associating 
local spin-polarisation axes with all lattice sites and the 
orientations {e^} vary very slowly on the time-scale of 
the electronic motions. [j| These 'local moment' degrees 
of freedom produce local magnetic fields on the lattice 
sites which affect the electronic motions and are self- 
consistently maintained by them. By taking appropriate 
ensemble averages over the orientational configurations, 
the magnetic properties of a system can be determined 
and, with the explicit inclusion of relativistic effects upon 
the electronic structure, the temperature dependence of 
the MCA can be obtained. 

Consider this DLM picture of a ferromagnetic metal 
magnetised along a direction h at a temperature T. The 
probability that the system's local moments are config- 
ured according to {e{\ is 

P^{{h}) = exp[-/?r!(")({e 4 })]/Z( ft ) (1) 
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where the partition function Z^"' = 
n,/rfe,exp[-/30( fi )({e l })]. nW({e<}) is the 'gen- 
eralised' electronic grand potential from SDFT Q and 
(3 = (ksT)^ 1 . The thermodynamic free energy, which 
accounts for the entropy associated with the orienta- 
tional fluctuations as well as creation of electron-hole 
pairs, is given by F*' = —ksTlogZ^. The role of a 
local moment hamiltonian, albeit a highly complicated 
one is played by By choosing a suitable reference 

'spin' hamiltonian Slnje,; ) = Y^i^ 1 ^ ■ e, and, using 
the Feynman Inequality [21|, a mean field theoretical 
estimate of the free energy is obtained, 

F (n) = + (i //3 )^ J pWfg,) lnpW (e*)^ (2) 

i 

where the probability of a moment pointing along ej on 
a site i is 



P {fl \e t ) = 



exp[—0h^'n ■ ei 



J exp[— fth^n ■ ei]de.i 
and the Weiss field at a site is given by 



(3) 



(4) 



where (• ■ ■}g i denotes a constrained statistical average 
with the magnetic moment on site i being fixed along ej . 

The magnetisation M = Mn is given by M = 
\x J P( n '(ei)n ■ iidii. fj, is the size of the local mo- 
ment on the site and is determined self-consistently. Q 
For materials where this DLM picture is suitable, the 
sizes of the local moments, \x, remain fairly constant 
so that even in the paramagnetic state where M = 0, 
the /x's are roughly the same as the magnetic moment 
per atom in the ferromagnetic state at T — OK. In a 
first-principles implemenation of such a DLM picture, 
the averaging over the orientational configurations of 
the local moments is performed using the KKR-CPA 
method adopted from the theory of random metallic al- 
loys. [(J [22] Using this methodology over the past 20 
years, the paramagnetic state, onset of magnetic order 
and transition temperatures of many systems have been 
successfully described. |2j| All applications to date, how- 
ever, have neglected relativistic effects and have been 
devoted to the paramagnetic state where the symme- 
try turns the calculation into a binary alloy-type one 
with half the moments oriented along a direction and the 
rest antiparallel. Once relativistic effects are included 
and/or the ferromagnetic state is considered, this sim- 
plicity is lost and the continuous probability distribution 
PW(ej)'s must be sampled for a fine mesh of angles and 
the averages with the probability distribution performed 
numerically. (Careful checks have to be made to en- 
sure that the sampling of P^ n \e.i) is sufficient - in our 
calculations some 25,000 values were used.) In the ferro- 
magnetic state, the magnetic anisotropy is given by the 
difference between the free energies, for different 



magnetisation directions, n, but the same magnetisation 
M. 

Once again our study of FePt starts with a self- 
consistent, scalar-relativistic calculation, this time for 
the paramagnetic (DLM) state. On the Fe sites a local 
moment of 2.97 hb is set up whilst no moment forms on 
the Pt sites. For the same lattice spacings (c = 0.385nm, 
c/a = 1), we found that, for the completely ferromag- 
netically ordered state of FePt at T = OK, the magneti- 
sation per Fe site is 2.93 [Ib and a small magnetisation 
of 0.29 [is is associated with the Pt sites. This sug- 
gests that the thermal effects on the magnetic properties 
should be well-described by the DLM picture. Using the 
self-consistent potentials and magnetic fields of the para- 
magnetic DLM state, we proceed by picking a series of 
values of M(T)/M(0) to set the probabilities, P {fl) {ei). 
By calculating (n (A) {e/})g 4 , using our relativistic DLM 
method, we find the Weiss field, h n , at each value and 
for a given n, and hence determine the temperature de- 
pendence of the magnetisation. The results are shown 
in Figure 1. Although the shortcomings of the mean 
field approach do not produce the spinwave T= behav- 
ior at low temperatures, the easy axis for the onset of 
magnetic order is obtained ( M 001 ) > ft/ 100 ) as T — > T c ) 
and it corresponds to that found at lower temperatures 
both experimentally and in all theoretical (T = OK) cal- 
culations, perpendicular to the Fe, Pt layer stacking. 
(A straightforward adaptation of this approach to other 
systems such as thin films in combination with T = OK 
calculations may be useful in understanding temperature 
induced spin reorientation transitions.) A Curie tem- 
perature of 935K is found, in fair agreement with the 
experimental value of 750K. [REif 
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FIG. 1: The magnetisation of ordered FePt as a function 
of temperature. The filled squares refer to a magnetisation 
along n — (0, 0, 1) whilst the open circles refer to h = (1, 0, 0). 
The onset of magnetic order occurs at 935K along the easy 
axis, (0, 0, 1). The full line shows the magnetisation from the 
mean field approximation to a classical Heisenberg model for 
comparison. The inset shows the region near the onset. 
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The free energy difference, F^°' '^ - i^ 1 * ' ), i.e. the 
magnetocrystalline anisotropy, K(T), is also calculated 
using the theory and leads to the key results of this let- 
ter shown in Figure 2. At T = OK, the MCA has a 
value -1.740 meV, close to the value, -1.696 meV, ob- 
tained by the earlier, separate calculation for this sit- 
uation for the completely ferromagnetic state. Figure 
2 shows the MCA, K{T), versus M(T)/M(0). Curves 
for the single-ion classical spin model anisotropy and 
(M (T) / M (Q)) 2 are also shown for comparison. Apart 
from near complete magnetisation, 0.9 < M(T)/M(0) < 
1, our results show a good approximate (M(T)/M(0)) 2 
behavior in good agreement with experiment. [T3 . Hq . 
This is in marked contrast to the model, which becomes 
proportional to (Af (T)/M(0)) 3 for the larger M(T)'s. 
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FIG. 2: The magnetic anisotropy of FePt as a function of 
magnetisation. The filled circles show the calculations from 
the ab-initio theory. The full line show Ko(-ff(jjy) 2 and the 
dashed line shows the single-ion model function Kn < . 32 [ e l^ T 
with K = -1.835 meV. 

Evidently Figure 2 shows that at low temperatures 
in the single ion model the MCA falls off much more 
quickly as the temperature is increased and the over- 
all magnetisation is reduced. Moreover our theory for 
an itinerant electron system does capture the behavior 
of the K versus magnetisation relation quantitatively. 
This theory assumes that there is a separation between 
fast and slow electronic degrees of freedom. A picture 
of 'local moments' emerges naturally but with a sub- 
tlety that their existence and behavior are determined 
by the fast electronic motions. We expect the MCA of 
the important magnetic materials Llo-CoPt and FePd 
to follow a similar variation with magnetisation since 
the local moments sustained on the Co and Fe sites in 
the paramagnetic DLM states (1.78 and 2.98 [Ib respec- 
tively) are comparable in size to magnetisation per site in 
the completely ferromagnetic states (1.91 and 2.96 [Ib)- 
Our DLM theory should therefore have good prospects 
in describing the variation of K with magnetisation for a 



range of metallic magnets like these. The success of the 
above relativistic DLM methodology in explaining the 
unexpected behavior of LIq- FePt suggests that further 
calculations for promising magnetic materials in bulk, 
thinfilms or in magnetic nanostructures may be valuable 
for the future modelling and exploitation of their mag- 
netic properties. 
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